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We present a detailed numerical study of a microscopic artificial swimmer realized recently by Dreyfus et 
al. in experiments [R. Dreyfus et al.. Nature 437, 862 (2005)]. It consists of an elastic filament composed of 
superparamagnetic particles that are linked together by DNA strands. Attached to a load particle, the resulting 
swimmer is actuated by an oscillating external magnetic field so that it performs a non-reciprocal motion in order 
to move forward. We model the superparamagnetic filament by a bead-spring configuration that resists bending 
like a rigid rod and whose beads experience friction with the surrounding fluid and hydrodynamic interactions 
with each other We show that, aside from finite-size effects, its dynamics is governed by the dimensionless 
sperm number, the magnitude of the magnetic field, and the angular amplitude of the field's oscillating direction. 
Then we study the mean velocity and the efficiency of the swimmer as a function of these parameters and the 
size of the load particle. In particular, we clarify that the real velocity of the swimmer is influenced by two 
main factors, namely the shape of the beating filament (determined by the sperm number and the magnetic-field 
strength) and the oscillation frequency. Furthermore, the load size influences the performance of the swimmer 
and has to be chosen as a compromise between the largest swimming velocity and the best efficiency. Finally, 
we demonstrate that the direction of the swimming velocity changes in a symmetry-breaking transition when 
the angular amplitude of the field's oscillating direction is increased, in agreement with experiments. 

PACS numbers; 87.19.St, 87.16.Ac,87.16.Qp 



I. INTRODUCTION 

Nature was very inventive to design mechanisms that mi- 
croorganisms such as bacteria and many eukaryotic cells use 
to propel themselves in a highly viscous environment, i. e., 
at very low Reynolds numbers [1]. Since they cannot rely on 
drifting by inertia, as we do when we swim in water, they im- 
mediately come to a halt when they stop with their beating 
motion. In this article we study in detail a microscopic artifi- 
cial swimmer that was constructed recently on the basis of 
a superparamagnetic elastic filament that mimics the so-called 
flagellum employed by many eukaryotic cells 1 1 ]. 

In 1977 Purcell pointed out in his famous article "Life at 
low Reynolds number" that microorganisms have to perform 
a non-reciprocal periodic motion to be able to move forward 
lit] (see also Non-reciprocal means that the time-reversed 
motion is not the same as the original one (for examples see 
llLP 0, H H [Ml). The reason Ues in the Stokes equations 
I IJI governing the fluid flow around the microorganisms for 
negligible inertia: they allow for a time-inverted flow pattern 
when all the external forces are inverted. 

Bacteria employ a marvelous rotary motor to crank a rela- 
tively stiff helical filament lfl2[ [Tsll . Spermatozoa as one ex- 
ample for eukaryotic cells move forward by creating bending 
waves that move along their elastic flagella from the head to 
the tail lfl4i[l5i[l6ll . While these waves are generated by a col- 
lective motion of internal molecular motors iflTlflsll . Dreyfus 
et al. use an external magnetic field to induce the beatingof 
a superparamagnetic filament attached to a red-blood cell H • 
The filaments are made from superparamagnetic colloidal par- 
ticles of micron size. A static external magnetic field induces 
dipoles in the colloids so that they form a chain. In the gaps 
between the charged colloids chemical linkers such as double- 
stranded DNA are attached to the particles and an elastic fil- 
ament resisting bending and stretching is formed 11^ I21I1 



(for similar systems see 112211 ). As demonstrated by the im- 
pressive experiments of Dreyfus et al, an oscillating external 
magnetic field now induces a non-reciprocal beating motion 
of the superparamagnetic filament that is able to move the at- 
tached red-blood cell forward. 



So far, the modelling of the dynamics of the superparamag- 
netic filament followed the elastohydrodynamics of an elastic 
rod ifisi I23I1 supplemented by a continuum version for the in- 
teraction of the magnetic-field induced dipoles H^l or a 
simpler description for the interaction with the magnetic field 
ll25ll . The authors are able to describe the dynamics of the 
filament ll24il25ll and the velocity curve of the artificial swim- 
mer 101 ■ Here we present a different description of the arti- 
ficial swimmer similar to the one used by Lagomarsino and 
Lowe for driven microfilaments ll26i l27ll . We take into ac- 
count the discrete nature of the superparamagnetic filament 
by modelling it as a seguence of beads, and, in contrast to the 
work of Roper et al. |24], we consider dipolar and hydrody- 
namic interactions between all the beads. We present a thor- 
ough investigation of the swimming velocity and efficiency as 
a function of the relevant parameters, study the influence of 
the size of the load attached to the filament, and demonstrate 
that the swimmer can take different directions via a symmetry- 
breaking transition depending on how the actuating magnetic 
field oscillates. 



In Sec. we present details for the modeling of the dy- 
namics of the superparamagnetic filament and show that aside 
from finite-size effects it is governed by a few relevant param- 
eters. Section Hill summarizes and discusses the results from 
our numerical study and Sec. |IV] contains our concluding re- 
marks. 
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FIG. 1 : Bead-spring model for the superparamagnetic filament con- 
sisting of N particles. Each bead at position r, experiences forces due 
to bending and stretching and due to interactions of the magnetic- 
field induced dipoles located on the beads. 



II. MODELLING THE SUPERPARAMAGNETIC 
FILAMENT 

We model the superparamagnetic filament by a bead-spring 
configuration, as illustrated in FigJT]that, in addition, resists 
bending like a worm-like chain luSll . Thus each bead in the 
filament experiences a force due to stretching, bending and 
dipolar interactions. While the chemical linkers between the 
beads are responsible for the resistance to bending and stretch- 
ing, we completely ignore their contribution to hydrodynamic 
friction. So the filament interacts with the fluid surrounding 
via flie hydrodynamic friction of the beads that incorporates 
their hydrodynamic interactions. In the following we first set 
up the discretized version of the bending and stretching free 
energies ll29ll and also formulate the dipolar interaction energy 
from which we then calculate the forces on a single bead. In 
a second step, the equations of motion for the single beads 
are established with the help of their mobilities on the Rotne- 
Prager level. Comments on numerical details and the scaling 
behavior conclude this section. 



A. Energies and Forces 

A deviation of the distance Z, ~ \ti\ of the beads from their 
equilibrium value Iq obeys Hooke's law and the total stretch- 
ing free energy is 



(1) 



1=2 



where denotess the total number of beads. Typically, we 
consider relatively stiff springs so the variation of Z,- is always 
smaller than 0. 1 /q. The stretching force acting on bead / obeys 

Ff = -Vr^H^ = -k{U - lo)t. + k{li+i - lo)ti+i , (2) 

where V^, is the nabia operator with respect to r,. 

The linkers connecting the beads convey some bending 
rigidity to the filament. The bending free energy of an elastic 
rod or a worm-like chain is given by llsoll 



R 1 fdt 



(3) 



where L is the total length of the filament, s the arclength 
along it, and t the unit tangent at location s. Refering the 
bending stiffness A = k^Tlp to the thermal energy k-^T , one 
obtains the persistence length Ip that gives the length scale on 
which the filament becomes flexible. Replacing dt/{ds) by 
(t,+ i — ii) /Iq, we arrive at the discretized version of the bend- 
ing energy for the superparamagnetic filament. 



A ^ 



where i, = We have introduced the factor 



fi 



1 for 2 < / <A^- 1 
for /= 1,N 
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to let the index ; in Eq. (|4|i still run from 1 to A^, i. e., over 
all beads. This facilitates the calculation of the bending force 
F^^ = —'VrjH^ acting on each bead ;; 
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To arrive at Ff we have used the relations V^,*/ = j: (1 — t,- 

£,) and Vr,t;+i = —7^(1 — ti+i ®ti+i) where 1 is the unit 
tensor and the symbol (g) means tensor product. 

A single particle made from material with magnetic suscep- 
tibility X and subjected to an external magnetic field B (more 
correctly the magnetic induction) aquires a magnetic dipole 
moment 



3a'o 



(7) 



where /jq = 471 • lO^^N/A^ is the permeability of free space 
and a is the particle's radius. In the filament, the dipole mo- 
ments Pi of the beads interact with each other with a total 
dipole-dipole interaction energy 
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where r,y = |ry — r,|, fij = {rj — ri)/rij and Y! means sum 
over all / 7^ j. When all the beads experience the same mag- 
netic field B = Bp, the total interaction energy becomes 



r, 4na^ , ^2 1 
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Here we ignore that the local magnetic field differs from the 
external B since all the induced dipoles contribute to the local 
field. Taking into account the field from the nearest neighbors 
of a dipole, one can show that this renormalizes the suscepti- 
bility X and gives it a small anisotropy with X\\ ~ X± ~ —0.25 
for a typical value of x ~ 1, where || and _L refer, respectively, 
to directions parallel and perpendicular to the local filament 
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axis 112411 . In the following we are interested in the basic fea- 
tures of the system and therefore work with the energy (|9]l to 
calculate the force = — Vr-,//^ that all magnetic dipoles 
excert on bead /; 



m ,=1 



rij)p+[l-5{p-fij)^]fi^^ 



(10) 



The filament is now driven by an external field whose direc- 
tion p{t) and magnitude B{t) depend on time. For example, if 
the direction is oscillating very slowly, the filament will fol- 
low the field like a rigid rod since the induced dipoles prefer to 
align themselves along a common line. However, if the oscil- 
lations become faster, the hydrodynamic friction of the beads 
with the surrounding fluid increases so that they are not fast 
enough to form a straight line and the filament is bending. 



To obtain the particle paths r,(f) and thus the dynamics of 
the filament, we numerically integrate Eqs. (fTTT i using the Eu- 
ler method Issll . More accurate schemes such as the Runge- 
Kutta method that would allow larger time steps during the 
integration and thus speed up the simulations are not useful 
since the maximum time step for the integration is governed 
by the relaxation dynamics of local bending and stretching de- 
formations. At the time step chosen to avoid numerical insta- 
bilities, the performance of the Euler method is comparable, 
e. g., to the fourth-order Runge-Kutta scheme. In concrete, the 
relative differences of the particle positions calculated with 
both methods are smaller than lO^"*. 



C. Reduced Equations of Motion 



B. Equations of Motion 

The filament is immersed in a viscous fluid such as water. 
On the length scale of microns and for times larger than the 
momentum relaxation time, inertia can be neglected so that 
the velocities v, of the beads are proportional to the forces 
acting on them fllll . Hence the beads follow the equations of 
motion 



Y^fiijFj with Fj 



-Ff- 



(11) 



where the forces F^ , F^, and Fj' are given in the last sec- 
tion. They all depend on the beadslocations r, and the dipolar 
forces also possess an exphcite time dependence. The impor- 
tant quantities in our treatment are the mobilities fiij. Roper 
et al. considered the filament as a continuous line whose fric- 
tion with the surrounding fluid is governed locally by two 
anisotropic friction coefficients for respective motions parallel 
and perpendicular to the local direction of the filament ll24ll . 
In our treatment, the anisotropic friction results from hydrody- 
namic interactions between neighboring beads. Moreover, hy- 
drodynamic interactions between more distant beads are also 
taken into account. Since induced flow fields are long ranged 
(they decay as 1 /r, where r is the distance from a moving 
bead), this effect cannot be neglected. In general, hydrody- 
namic interactions constitute a complicated many-body prob- 
lem ifTHl . however, their leading contribution is given by two- 
particle interactions. Moreover, if the particles are not too 
close to each other so that lubrication becomes important, the 
Rotne-Prager approximation can be employed ifTu [3lll . We 
use it in a version for spheres with different radii a, and cij 
isill . Whereas all spheres in the filament have the same radius 
at = a {\ < i < N), we also attach a larger sphere with ra- 
dius flo to the filament to mimic its load. In the Rotne-Prager 
approximation, the self mobility is simply 



flu = ^ol with l/(67irifl,) , 



(12) 



and the cross mobilities read 
1 
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Our modelling of the filament comprises several parame- 
ters. To identify the essential parameters that govern the dy- 
namics of the filament, we rescale the dynamic equations ap- 
propriately so that only reduced variables appear For exam- 
ple, all lengths will be referred to the equilibrium length of the 
filament that we approximate by L « NIq. 

To arrive at reduced equations of motion, we first rescale 
the energies from section UTaI with the help of characteristic 
quantities: 
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where the reduced energies read: 
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(15) 
(16) 
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The prefactor kL^/N of in Eqs. ( fT4] i is essentially the 
stretching modulus of an elastic rod and averages the 
square of the strain variable over the whole rod. The bending 
free energy of the filament with a uniform curvature L/ \f7. 
amounts \a AjL and the prefactor of ff^ gives the order of 
magnitude of the energy of interacting dipoles. To compare 
these characteristic energies to each other, we introduce the 
reduced magnetic field Bj ll36ll and the parameter ks. 



■4\ 1/2 



4na 



2N-* 



A/L 



kL-/N 
A/L 



2,3 



B (18) 



(19) 



In the following, we will characterize the strength of the mag- 
netic field by B,, the parameter kg will basically be kept con- 
stant. 

Whereas the scaling factors in Eqs. (fT4l i determine the mag- 
nitude of the energies, the reduced energies of Eqs. (fT5]l-(fT7l) 
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distinguish between different types of stretching and bending 
deformations and arrangements of dipoles along the filament. 
For sufficiently large L « NIq, they tend to constant values 
H^,H^, and h2- We checked how finite size effects from the 
ends of the filaments influence these constant values. Bending 
the filament to a semicircle, the deviation from is around 
1 % when = 20, a typical number of beads in the super- 
paramagnetic filament both in experiments and in our simula- 
tions. On the other hand, for the same number of beads 
for dipoles aligned along a straight line deviates from h2 by 
around 10% due to the long range nature of the dipolar inter- 
action. 

To formulate the reduced equations of motion, we rescale 
length, time, and velocity according to 



, . df \ dr 

f = r L , f = (Of , — = , 

dt &L dt 



(20) 



where CO is the frequency of the actuating magnetic field. The 
reduced forces and mobilities follow from 



(21) 



Note that 6Kr[aL/lQ is the total friction coefficient of the fila- 
ment without taking into account hydrodynamic interactions. 
Applying the rescaled quantities to Eqs. ( fTTT ) and using Eqs. 
(O, ( fTSl ), and (T% . one arrives at the reduced equation of 
motion for bead /: 



dfi 
dt 



with 



Vrj{H" + k,Hy2+B;H 



(22) 



(23) 



The important parameter Sp in Eqs. (l22l i. introduced first by 
Wiggins and Goldstein ll23l IjtIi and termed sperm number by 
Lowe ll26ll . is defined via 
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It compares the frictional to the bending forces and completely 
determines the dynamics of an elastic filament in a viscous 
environment in the so-called resistive-force theory of slender 
bodies ifssll . where the friction with the surounding fluid is de- 
scribed by two local friction coefficients per length, y|| and yx, 
for respective motions parallel and perpendicular to the local 
axis of the filament. Our definition of Sp agrees with the one 
of Wiggins and Goldstein ll23l IstIi when we make the reason- 
able identification yx = 67tr|fl /Iq. An immediate interpretation 
of Sp —L/li^ is given with the help of the elastohydrodynamic 
penetration length //, ||23[ IstIi . Consider a sufficiently long fil- 
ament (L 3> //,) whose one end undergoes an oscillation with 
frequency co. Then Z/, is the length on which the oscillation 
penetrates into the filament. On the other hand, if L <C //,, the 
filament oscillates as a whole like a rigid rod. 

Other important quantities are the strength of the actuat- 
ing magnetic field quantified by the parameter Bs and its time 



protocol, which we will concretize in section|III] The reduced 
stretching constant is always chosen such that variations in 
the length of the filament are smaller than 10%. 

We already mentioned that the dynamics of a continuous 
elastic filament described within the resistive-force theory is 
completely determined by the sperm number Sp, which espe- 
cially incorporates the influence from the filament length. In 
the superparamagnetic filament simulated by us through a dis- 
cretized model, finite-size effects are present, as already dis- 
cussed for the bending and dipolar energies. Hydrodynamic 
interactions are of longer range than dipolar interactions and 
we therefore discuss shortly how variations in the length of 
the filament or in the number of beads influence the dynamics 
for constant Sp, Bs, and ks- We keep the strength of the mag- 
netic field constant but let its direction oscillate in the yz plane 
around the z axis with an angular amplitude of 50°. We solve 
the equations of motion ( fTTT i for real parameters and change 
the length of the filament for constant equilibrium-bead dis- 
tant lo by choosing the particle numbers = 20, 40, and 80. 
The reduced parameters are kept fixed at the values Sp = 7.7, 
Bs = 6.6, and ks = 4500 by changing, respectively, the fre- 
quency CO of the oscillating magnetic field, its strength B, and 
the real spring constant k. Figure |2^) shows snapshots of the 
different filaments in the yz plane at the same moments within 
the period of the oscillating field. Hydrodynamic interactions 
are completely switched off. For = 40 and 80 the con- 
figuration of the filaments are basically identical whereas for 
N — 20 one realizes small differences due to the finite-size ef- 
fects from the dipolar interactions as discussed above. In Fig. 
I2J3) the configurations are shown at the same moments within 
the period as in Fig. [2^) but now with hydrodynamic interac- 
tions switched on. They clearly have a pronounced effect on 
the oscillating shape of the filament. Since the beads move 
collectively through the fluid, hydrodynamic interactions re- 
duce, roughly speaking the local friction with the surround- 
ing fluid. The filament can more easily follow the oscillating 
magnetic field direction and, as a result, it is less bent. The 
finite-since effect is more pronounced compared to Fig. |2^) 
but does not change the configurations significantly. So when 
in the following all our simulations are done with = 20 to 
reduce computing time, we are sure that the dynamics of the 
filament is mainly described by the sperm number Sp, the re- 
duced strength Bs of the magnetic field, and its time protocol. 



III. STUDY OF THE ARTIFICIAL SWIMMER 

The actuated superparamagnetic filament in Fig.|2]does not 
move on average in one direction since it performs a recipro- 
cal motion. This symmetry is broken when a load is attached 
to one end of the filament as demonstrated by the wonderful 
experiments of Dreyfus et al. where the load is a red-blood 
cell In the following we study such an artificial swim- 
mer in detail. As illustrated in Fig. [3j we actuate the fila- 
ment with a magnetic field B{t) whose strength is constant 
but whose direction oscillates about the z axis with an angle 
(p(f) = (pniaxsin(cof) (this time protocol differs from the one 
used in Ref. As a result, the swimmer will move with 
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FIG. 2: Configurations of the superparamagnetic filament at the same 
moments within one period of the oscillating magnetic field. The 
reduced parameters Sp,Bs, and are kept constant while the number 
of beads is changed, a) without hydrodynamic interactions, b) with 
hydrodynamic interactions. The numbers at the snapshots indicate 
the pair of variables (p, cof / (27l) for the time protocol of the oscillating 
magnetic field, as illustrated in Fig.[3l (1) 50°, 0.25; (2) 0°,0.5; (3) 
-35°, 0.88. 



an average velocity v along the z axis. However, in contrast 
to spermatozoa, where the head is pushed forward by damped 
waves travelling from the head to the tail lfl4i[l5ll . the super- 
paramagnetic filament drags the passive load behind itself by 
performing a sort of paddle motion with its free end as indi- 
cated in Fig. |3] Note that the direction can be reversed for 
certain parameters when the load particle also becomes super- 
paramagnetic, but we will not study this possibihty further. 
In our approach, hydrodynamic interactions between the con- 
stituent beads are essential for the artificial swimmer to move 
forward. In particular, they produce local friction coefficients 
Y|| and Yx for respective motions parallel and perpendicular 
to the filament's local axis that are not equal, Y|| ^ Y±- This 
is a necessary prerequisit for swimming ||27I1 . To implement 
hydrodynamic interactions, we use the Rotne-Prager approx- 
imation for the mobilities as introduced in Eqs. (fT2] i and (flJT l 
in the form that also allows for hydrodynamic interactions be- 
tween the filament's beads and the larger load particle. We 
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FIG. 3: For the artificial swimmer, a larger load particle is attached 
to the filament. It is actuated by a magnetic field whose direction 
is oscillating about the z axis with (p(f) = (pniax sin(co/). The two 
configurations of the swimmer are schematic drawings for positive 
and negative (p(/). 



simulate the artficial swimmer by solving the equations of 
motion ( fTTI ) numerically for two sets of parameters summa- 
rized in Appendix |A1 The first set is close to the realization 
of the artificial swimmer in Ref. Since we already now 
that the filament itself is governed by the dimensionless pa- 
rameters Sp, Bs, and ks, we present our results as functions 
of the sperm number Sp and the reduced field strength Bs- In 
addition, we also show that the angular amplitude (pmax of the 
magnetic-field oscillations and the size oq of the load particle 
influence the dynamics of the swimmer. 

We discuss the performance of the swimmer by studying in 
detail its average speed -0 and its efficiency to move a load. 
The speed v is calculated by averaging the velocity vq of the 
load particle over one actuation cycle: 



vo{t)dt , 



(25) 



where T = 2n/(>i. In order to obtain a reliable v from Eq. ( |25] ), 
it is important that the swimmer has reached a steady state at 
time T so that the swimming does not depend any more on 
initial conditions. For small frequencies {Sp « 1), this state 
is reached rather quickly after one actuation cycle but needs 
several cycles when Sp is increased. Secondly, we define the 
efficiency of the swimmer to transport a load by comparing 
the energy dissipated by the load, when moved uniformly with 
velocity v, to the total energy dissipated by the swimmer: 



67tr|flov 



n 



(26) 



where the bar in the denominator means average over one ac- 
tuation cycle. The efficiency t, indicates how much energy 
from the total energy used to actuate the swimmer is employed 
to move the load particle forward with velocity v. 

In Fig. m we plot the mean velocity v ~ \v\ and the effi- 
ciency ^ as a function of the sperm number Sp. For the sim- 
ulations, the realistic parameter set I was used. The curve for 
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\'l (Leo) in the upper graph o f Fig . |4]follows the same behav- 
ior as observed in Ref. ll26l l27ll . where the elastic filament 
is driven by an oscillating torque acting on one of its ends. 
For small sperm numbers around 5,, = 3, the reduced velocity 
is small since the superparamagnetic filament behaves nearly 
Uke a rigid rod, as illustrated by the inset on the lower left. 
It shows snapshots of the filament for 5,, = 3. As already 
mentioned, the oscillating motion of a rigid rod is reciprocal 
and therefore cannot produce a directed motion of the swim- 
mer. Increasing the sperm number, e. g., by increasing the 
frequency CO, speeds up the artificial swimmer. At the maxi- 
mum value of \'l (Leo) at around Sp = 6, the increased friction 
with the surrounding fluid is able to bend the whole filament 
(see upper inset) which obviously promotes a high swimming 
velocity. Further increase in Sp leads to a decrease in v/ (Leo); 
due to the strong friction with the fluid the whole filament can- 
not follow the magnetic field and only a small wiggling of its 
free end remains. The efficiency as a function of Sp shows 
a similar behavior as v/(Lco) that is not completely surpris- 
ing: oscillating a rigid rod (small Sp) or fast wiggling of the 
filament (large Sp) dissipates energy but does not produce an 
effective motion. So one expects a maximum of ^ close to the 
maximum of v/ (Leo) since ^ is determined by v^. Note that 
the small absolute efficiency (^max = 1 -58 • 10^^) is due to the 
fact that a large amount of energy is dissipated by the motion 
of the filament. The shape of the velocity curve changes when 
absolut velocities are plotted, as demonstrated in the lower 
graph of Fig.m At Sp = 3, the absolute velocity is nearly zero 
(due to the small frequency) and the maximum is shifted to a 
larger value around Sp = 7.5. Interestingly, the absolute veloc- 
ities of the oscillating filaments at Sp ^6 and 1 2 do not differ 
so much, as first implied by the upper graph. This suggests 
that the absolute swimming velocity depends on two factors: 
(1) the shape of the oscillating filament, where bending the 
whole filament favors large velocities, and (2) the oscillation 
frequency. The latter point leads to the very slow decrease of v 
as a function of Sp in the lower graph. A comparison with the 
narrow maximum of ^, however, shows that a lot of energy is 
dissipated at large Sp in the small wiggling motion of the fil- 
ament. So operating the artificial swimmer at around Sp ~1 
between the two close maxima ensures highest swimming ve- 
locities with very efficient energy consumption. 

In Fig. |5]we present reduced velocities and efficiencies as 
a function of Sp for several reduced field strengths B,, which 
we simulated with parameter set II. We observe that with in- 
creasing Bs the maxima of both the velocity and the efficiency 
curves move to larger sperm numbers. We understand this 
since larger magnetic fields mean stronger alignment of the 
dipoles and, therefore, larger resistance to bending. Both, 
maximium velocity and efficiency, only exhibit a small de- 
pendence on Bs- After a slight increase of the maximum ve- 
locity, it slightly decreases to a constant value (not shown in 
the graph). Note that the reduced velocities are larger by a 
factor of 1.2 compared to Fig.|4]that we attribute to the larger 
angular amplitude (pmax = 57° compared to cpmax = 45° used 
in Fig.|4] Finally, Fig.|6]summarizes the dependence of the ab- 
solut swimming velocity on sperm number and reduced field 
calculated with parameter set II. For increasing field strength 
at constant Sp, the graph shows a strong increase of the swim- 




sperm number Sp 

FIG. 4: Swimming velocity v and efficiency ^ in units of ^max — 
1.58 • 10^^ as a function of sperm number Sp for reduced magnetic- 
field strength Bs = 5.76. Upper graph: reduced velocity v/{Lco), 
lower graph: absolut velocity v in units of 5.56- lO^V/s. 

The insets show several snapshots of the filament's configuration for 
Sp = 3,6 and 12, respectively, indicated by the dots. The numbers 
at the snapshots indicate the pair of variables (p, cof/ (2jl) for the time 
protocol of the oscillating magnetic field, as illustrated in Fig.|3] (1) 
20°, 0.07; (2) 40°, 0.17; (3) -20°, 0.57; (4) -40°, 0.67. Parameter 
set I was used in simulations. 



ming velocity. At small Sp a subsequent saturation at a con- 
stant value is visible. Note that the maximum velocity for con- 
stant Bs is shifted to higher sperm numbers when Bs increases. 
This explains why the the maximum reduced velocity v / (Leo) 
in Fig.|5]does not exhibit a strong dependence on By. 

The swimming velocity v clearly depends on the size or the 
radius aq of the load particle. In the upper graph of Fig. |7] 
we plot its absolute value as a function of sperm number Sp 
and flo for a constant field strength Bs = 5.76 calculated with 
parameter set I. There is a pronounced maximum at Sp = 8 
and flo ~ 3fl (indicated by a filled circle), where a is the ra- 
dius of the superparamagnetic beads in the filament. The 
velocity v decreases for large flo since the load becomes too 
heavy to be efficiently moved by the oscillating filament. On 
the other hand, when flo approaches a, we also expect small 
swimming velocities since the asymmetry of the swimmer be- 
comes small. However, even at flo = fl the swimming ve- 
locity is not zero and therefore the swimmer still performs 
a non-reciprocal motion. The reason is that the load particle 
is not superparamagnetic and, therefore, a small asymmetry 
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sperm number 5, 



FIG. 5: Reduced swimming velocity v/(L(£)) and efficiency ^ in units 
of ^max = 1.87 • 10^-' as a function of sperm number Sp for several 
reduced magnetic-field strengths: (1) 6.5 = 4.2, (2) Sj = 6.07, and 
(3) Bs = 7.93. Parameter set II was used in simulations. 



velocityv 
(in units of v^ax) 




sperm number 5 



FIG. 6: Absolute swimming velocity v in units of the maximal value 
Vmax as a function of sperm number 5,, and reduced magnetic-field 
strength Bs. Parameter set II was used in simulations. 

remains. The plot for the reduced velocity v/ (Leo) looks sim- 
ilar to the lower graph of Fig. [Tj the absolute maximum is 
not that pronounced and the maxium for constant oq moves 
to smaller sperm numbers when oq is increased. The lower 
graph shows how the efficiency t, behaves as a function of Sp 
and flQ. The absolute maximum, indicated by an open circle, 
is at Sp = 6.6 and ao w 5fl. For comparison the location of 
the maximum of the swimming velocity is again shown by the 
filled circle. Since it is relatively sharp, one has to choose 
a compromise for the performance of the swimmer between 
the largest swimming velocity and the best efficiency. The 
location of the maximum velocity certainly depends on the 
magnetic-field strength. We expect that it is shifted to larger 
radii ao when Bs increases. However, we have not studied this 
dependence in detail. 

So far, the angular amplitude (pmax of the oscillating field 
was always chosen sufficiently small so that the swimmer was 



velocity v 
(in units of Vmax ) 

1 r 




efficiency t, 




FIG. 7: Swimming velocity v (in units of 7.31- 10"^m/s) and 

efficiency t, (in units of ^max = 1,54- 10^^) as a function of sperm 
number 5;, and load size qq in units of a. The filled and open circles 
indicate, respectively, the absolute maxima of v and ^. Bs = 5.76 and 
parameter set I was used in simulations. 



moving along the z axis. For large (pmax, however, the mean 
velocity v assumes a non-zero angle \\t with the z axis. The up- 
per graph of Fig. [8] summarizes our results. For Sp larger than 
7 and angular amplitudes beyond (pmax = 70°, the swimming 
angle \\t jumps from 0° to 90° and the swimmer thus moves 
perpendicular to the z axis. Such an abrupt transition was also 
observed in the experiments by Dreyfus et al. Ii39il . It is asso- 
ciated with a broken symmetry since the swimmer could also 
move into the opposite direction with \|/ = —90° depending on 
the initial condition. In the beginning the filament and mag- 
netic field point along the z axis. When the magnetic-field 
direction turns towards the upper-half space {y > 0), the fila- 
ment follows. For large oscillation frequencies co (large Sp), 
it cannot, however, follow the field direction into the lower- 
half space {y < 0) and \|/ = 90° results, as pictured in the third 
example of Fig. |9l For decreasing sperm number, this sharp 
transition becomes smoother involving swimming directions 
inclined relativ to the z axis. Finally for Sp < 3, the swim- 
mer always swims along the z axis by following the oscillat- 
ing field as a nearly rigid rod with only small bending. The 
simulations of the swimmer's movements for \(/ 7^ are very 
time consuming since the swimmer needs some time to reach 
a steady state. Several snapshots of its configurations for dif- 
ferent parameters are summarized in Fig.|9] their locations in 
the graphs of Fig. [8] are indicated by filled black circles. The 
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iPmax [deg] 




FIG. 8: Swimming velocity v (in units of 5.57- 10" V/s) and 

the swimming angle \|f (measured relative to the z axis) as a function 
of sperm number Sp and angular amplitude (pmax of the oscillating 
field. Bs = 3.81 and parameter set II was used in simulations. 



temporal evolution of the first and second configurations in 
Fig. |9] are the dynamic analog of the static hairpin structures 
that form when the aligning magnetic field is suddenly turned 
around by 90° lfl9|] . Finally, the lower graph of Fig.[8]plots the 
swimmming velocity v as a function of Sp and (pmax- Note that 
the Sp and (pmax axes are differently oriented compared to the 
upper graph. When the angular amplitude (pmax increases from 
0°, we observe that v also increases (only partially shown in 
the graph). At Sp = 7 and around cpmax = 70°, a sharp decrease 
in v" coincides with the abrupt change of the swimming direc- 
tion from x|/ = 0° to 90° in the upper graph. For decreasing 
Sp, the sharp drop in v becomes smaller due to the occurence 
of inclined swimming directions. 



IV. CONCLUSIONS 

In this article, we have presented a detailed numerical study 
of a microscopic artificial swimmer realized recently by Drey- 
fus et al. in experiments jiZ]. The elastic superparamagnetic 
filament is modeled by a bead-spring configuration that also 
resists bending via a disrete version of the worm-like chain. 
Friction with the surrounding fluid is described by the sin- 
gle beads that also experience hydrodynamic interactions with 




FIG. 9: Snapshots of the swimmer's configurations for different pa- 
rameters giving rise to different swimming directions. Their loca- 
tions in the graphs of Fig.[8]are indicated by filled black circles. The 
numbers at the snapshots indicate the pair of variables (p, CO// (27t) for 
the time protocol of the oscillating magnetic field. 5^ = 4.65: (1) 
0°,0.00; (2) 95°, 0.20; (3) 58°, 0.40; (4) -58°, 0.60; (5) -95°, 0.80. 
Sp = 5.34: (1) 0°,0.00; (2) 100°, 0.25; (3) 0°,0.50; (4) -100°, 0.75. 
Sp = 5.98: (1) 41°, 0.06; (2) 118°, 0.22; (3) 79°, 0.39; (4) -37°, 0.55; 
(5) -117°, 0.72; (6) -82°, 0.88. 



each other The swimmer composed of the filament and an at- 
tached load is actuated by a magnetic field whose direction 
oscillates. 

We show that the superparamagnetic filament, aside from 
finite-size effects, can be described by the dimensionless 
sperm number, the magnitude of the magnetic field, and the 
angular amplitude of the field's oscillating direction. We then 
study the mean velocity and the efficiency of the swimmer, 
which we define appropriately, as a function of these parame- 
ters and the size of the load particle. In particular, we clarify 
that the real velocity of the swimmer depends on two main 
factors namely the shape of the beating filament and the os- 
cillation frequency. For given magnetic-field strength an opti- 
mum sperm number (or oscillation frequency) can be chosen 
such that mean velocity and efficiency are close to their maxi- 
mum values. Whereas the maximum rescaled velocity v/(La)) 
as a function of the sperm number only exhibits a weak depen- 
dence on the magnetic-field strength B^, the real maximum ve- 
locity V strongly increases with Bs since its location is shifted 
to larger Sp °^ CO^/^. A study of the influence of the load size 
for a particular field strength reveals, the optimum load has 
to be chosen as a compromise between the largest swimming 
velocity and the best efficiency. For increasing angular am- 
plitude of the field's oscillating direction, the direction of the 
swimming velocity changes in a symmetry-breaking transi- 
tion that is sharp for large sperm numbers, becomes smoother 
for decreasing Sp, and ultimately vanishes around Sp = 3. Ac- 
cordingly, the jump in the swimming angle relative to the sym- 
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1.0 




FIG. 10: Swimniing velocity v and efficiency ^ in units of ^max — 
1.8 • 10^"^ as a function of sperm number Sp. Parameters are the 
same as the one used for the red experimental data points of Fig. 4 
in Ref. 10]. The experimental points are included in this figure with 
the symbol +. 

metry axis decreases from 90° to 0°. 

We also applied the time protocal of Ref. 101 (i.e., a constant 
z component and an oscillating y component of the magnetic 
field) to actuate the one-armed swimmer For the parameters 
of the red experimental data points of Fig. 4 in Ref. 101, we 
find nearly quantitative agreement, as illustrated in Fig. [TO] 
To achieve this, we had to rescale the actuating magnetic field 
by a factor of 2.53 to account for the larger distance of the 
beads in our modelling and therefore to compensate for the 
weaker dipole interaction compared to the swimmer in Ref. 
|l2tl. This makes sense since the strength of the dipole inter- 
action of neighboring dipoles is the important parameter for 
the actuation of the swimmer. Deviations between our simu- 
lations and the experimental results might be due to the fact 
that we use a spherical load particle compared to the oblate 
shape of the red blood cell used in Ref. 101 and that we neglect 
corrections to the actuating external field due to the induced 
dipole fields. 

In Ref. ll4(ll . the authors envisage micro machines moving 
in blood vessels and doing necessary repair work. The one- 
armed swimmer offers an interesting possibility to propel such 
micro machines. Our numerical studies demonstrate that theo- 
retical modelling helps to elucidate the basic features but also 



to optimize the performance of such micro machines. 
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TABLE I: Parameter sets I and II used in the numerical studies. 



parameter set I set II 

N 20 20 

a\pim] 0.5 2.25 

ao 5a 8a 

/q 3a 3a 

X 0.993 1.704 

Ti[Ns/m-] 10"^ IQ-^ 

k[N/m] 1.5-10"^ 1.0 

A[Nm] 4.5-10^22 6.75-10"'^ 

(£)[2%/s] 1...2285 100 ...9.4- 10"* 

^ Sp 2.903 . . . 20.074 3.636 ... 20.134 

S[T] 0.07 0.01 . . . 0.49 

^Sj 5.76 0.233 ...11.433 

9max[ °] 45 57 

At[s] IQ-^ IQ-S 

n, 15 15 



APPENDIX A: PARAMETER SETS 

In table |l] we summarize the parameters in set I and II 
used for our numerical studies. The oscillation frequency 
CO = 2k /T and the magnetic-field strength B were varied to 
study the respective ranges of sperm number Sp and reduced 
magnetic field B^. The spring constant k, the time step At for 
the Euler integration, and the number n.y of simulation cycles 
were adjusted as necessary. Typical values are shown. 



[1] D. Bray, Cell Movements: From Molecules to Motiliy, 2nd ed., 

Garland Publishing, New Yorlc (200 1 ) . 
[2] R. Dreyfus, J. Baudry, M. L. Roper, M. Fermigier, H. A. Stone, 

and J. Bibette, Nature 437, 862 (2005). 
[3] E. M. Purcell, Am. J. Phys. 45, 3 (1977). 
[4] H. C. Berg, PNAS 93, 14225 (1996). 

[5] A. Shapere and F Wilczek, Phys. Rev. Lett. 58, 2051 (1987). 
[6] H. A. Stone and A. D. T. Samuel, Phys. Rev. Lett. 77, 4102 
(1996). 

[7] L. E. Becker, S. A. Koehler, and H. A. Stone, J. Fluid Mech. 
490, 15 (2003). 



[8] J. E. Avron, O. Gat, and O. Kenneth, Phys. Rev Lett. 93, 
186001 (2004). 

[9] A. Najafi and R. Golestanian, Phys. Rev. E 69, 062901 (2004). 
[10] R. Dreyfus, J. Baudry, and H. A. Stone, Eur. Phys. J. B 47, 161 
(2005). 

[11] J.K.G. Dhont, An Introduction to Dynamics of Colloids (Else- 
vier, Amsterdam, 1996). 

[12] H. C. Berg and R. A. Anderson, Nature 245, 382 (1973). 

[13] H. C. Berg, E. coli in Motion, Springer Verlag, New York 
(2004). 

[14] G. Taylor, Proc. R. Soc. A 209, 447 (1951). 



10 



[15] J. Gray and G. J. Hancock, J. Exp. Biol. 32, 802 (1955); P. Satir, 
J. Cell Biol. 39, 77 (1968). 

[16] R. W. Linck, Cilia and Flagella in Encyclopedia of Life Sci- 
ences, Nature Publishing Group (2001), www.els.net. 

[17] S. Camalet, F. Julicher, and J. Prost, Phys. Rev. Lett. 82, 1590 
(1999). 

[18] S. Camalet and E Juliclier, New J. Phys. 2, 24 (2000). 

[19] C. Goubault, P. Jop, M. Eermigier, J. Baudry, E. Bertrand, and 

J. Bibette, Phys. Rev. Lett. 91, 260802 (2003). 
[20] L. Cohen-Tannoudji, E. Bertrand, L. Bressy, C. Goubault, J. 

Baudry, J. Klein, J.-E. Joanny, and J. Bibette, Phys. Rev. Lett. 

94, 038301 (2005). 
[21] A. Koenig, P. Hebraud, C. Gosse, R. Dreyfus, J. Baudry, E. 

Bertrand, and J. Bibette, Phys. Rev. Lett. 95, 128301 (2005). 
[22] S. L. Biswal and A. P Cast, Phys. Rev. E 68, 021402 (2003); 

Phys. Rev. E 69, 041406 (2004). 
[23] C. H. Wiggins and R. E. Goldstein, Phys. Rev. Lett. 80, 3879 

(1998). 

[24] M. L. Roper, R. Dreyfus, J. Baudry, M. Eermigier, J. Bibette, 

and H. A. Stone, to be published in J. Eluid Mech. (2006). 
[25] A. Cebers, J. Phys.: Condens. Matter 15, S1335 (2003); Curr. 

Opin. Coll. Interface Science 10, 167 (2005); Magnetohydro- 

dynamics 41, 63 (2005). 
[26] C. P Lowe, Phil. Trans. R. Soc. Lond. B 358, 1543 (2003). 
[27] M. C. Lagomarsino, E. Cap uani, C. P. Lowe, J. Theor. Biol. 

222, 215 (2003). 

[28] O. Kratky and G. Porod, Reel. Trav. Chim. Pays-Bas 68, 1106 
(1949). 



[29] Since the filament is actuated by an external field that only os- 
cillates in a plane, twisting of the filament does not occur and, 
therefore, does not have to be taken into account when formu- 
lating the free energy of the system. 

[30] L. D. Landau and E. M. Lifschitz, Lehrbuch der Theoretischen 
Physik Band VII: Elastizitdtstheorie (Akademie Verlag, Berlin, 
1991) 

[31] In our modelling, we use a center-to-center distance 3a of the 
particles in the filament (where a is the particle radius). In this 
case, the Rotne-Prager approximation deviates from the full 
multipole results of Ref. 1331 by at most 5%. We checked this 
by using the program Hydrolib (see Ref. (33ll ). 

[32] B. Cichocki, B. U. Eelderhof, K. Hinsen, E. Wajnryb, and J. 
Blawdziewicz, J. Chem. Phys. 100, 3780 (1994). 

[33] K. Hinsen, Comput. Phys. Commun. 88, 327 (1995). 

[34] D. J. Jeffrey and Y. Onishi, J. Eluid Mech. 139, 261 (1983). 

[35] W. H. Press, S. A. Teukolsky, W. T. Vetterling, and B. P Elan- 
nery. Numerical Recipes in C: The Art of Scientific Computing, 
Cambridge University Press, Cambridge (1992). 

[36] 5? is p roportional to the magnetoelastic number used in Refs. 

mm- 

[37] C. H. Wiggins, D. Riveline, A. Ott, and R. E. Goldstein, Bio- 

phys. J. 74, 1043 (1998). 
[38] C. Brennen and H. Winet, Ann. Rev. EluidMech. 9, 339 (1977). 
[39] R. Dreyfus, private communication (2005). 
[40] K. Ishiyama, M. Sendoh, A. Yamazaki, M. Inoue, K. I. Aral, 

IEEE Trans. Mag. 37, 2868 (2001). 



